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ABSTRACT 

Long-duration gamma-ray bursts (GRBs) are powerful tracers of star-forming 
galaxies. We have defined a homogeneous subsample of 69 Swift GRB-selected 
galaxies spanning a very wide redshift range. Special attention has been devoted 
to making the sample optically unbiased through simple and well-defined selec- 
tion criteria based on the high-energy properties of the bursts and their positions 
on the sky. Thanks to our extensive follow-up observations, this sample has now 
achieved a comparatively high degree of redshift completeness, and thus provides 
a legacy sample, useful for statistical studies of GRBs and their host galaxies. 
In this paper we present the survey design and summarize the results of our ob- 
serving program conducted at the ESO Very Large Telescope (VLT) aimed at 
obtaining the most basic properties of galaxies in this sample, including a catalog 
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of R and K s magnitudes and redshifts. We detect the host galaxies for 80% of the 
GRBs in the sample, although only 42% have i^-band detections, which confirms 
that GRB-selected host galaxies are generally blue. The sample is not uniformly 
blue, however, with two extremely red objects detected. Moreover, galaxies host- 
ing GRBs with no optical/NIR afterglows, whose identification therefore relies 
on X-ray localizations, are significantly brighter and redder than those with an 
optical/NIR afterglow. This supports a scenario where GRBs occurring in more 
massive and dusty galaxies frequently suffer high optical obscuration. Our spec- 
troscopic campaign has resulted in 77% now having redshift measurements, with 
a median redshift of 2.14 ±0.18. TOUGH alone includes 17 detected z > 2 Swift 
GRB host galaxies suitable for individual and statistical studies — a substantial 
increase over previous samples. Seven hosts have detections of the Lya emis- 
sion line and we can exclude an early indication that Lya emission is ubiquitous 
among GRB hosts, but confirm that Lya is stronger in GRB-selected galaxies 
than in flux-limited samples of Lyman break galaxies. 

Subject headings: cosmology: observations - galaxies: fundamental parameters 
- galaxies: high-redshift - gamma rays: bursts - surveys 



Introduction 



Following the discovery o f long-duration gamma-ray burst (GRB) optical afterglows 



[e.g., 



van Paradijs et al.l 120001 ) , evidence rapidly a ccumulated pointing to their origin in 



1998; Bloom et al 



2002; Stanek et al. 2003 



the deaths of ma ssive stars (IGalama et al 
Hjorth et al.l 120031 ) . This soon led to their being touted as a potentially powerful means 
of discovering and reckoning star formation acti yity, which could b ypass many of the biases 
that hamper other star formation measures (e.g.. lWijers et al.lll998l ). The detection of GRBs 
is at once both independent of the luminosity of their host galaxy and, at high energies, un- 
affected by dust obscuration. They are also extremely luminous at all wavelengths, so that 
they can be used to perform transmission spectroscopy of t heir host galaxies, even where 



such hosts are too faint to be detected photometrically (e.g.. lThone et al.l I201CH ) . This also 



permits them to be found across the universe from very low to very high redshift: indeed 



1 Based on observations collected at the European Southern Observatory, Paranal, Chile (ESO Large Pro- 
gramme 177.A-0591). Reduced data and catalogs are made available at http: //dark-cosmology. dk/TOUGH 
and http : / / archive . eso . org. 
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GRB 090423 h olds the record at z = 8.2 for the most distant spectroscopically confirmed 



object known ( ITanvir et al.ll2009i ; ISalvaterra et al.ll2009l ) 



The history of baryonic structure in the universe to a large extent revolves around the 
co-evolution of stars, galaxies, metals, and dust. In recent years, a range of methods has been 
applied to the identification and study of different populations of high-redshift galaxies, but 
all are variously affected, to a greater or lesser extent, by the problems of dust extinction, 
flux-limited samples, source confusion, an d incomplete redsh ift determinations. Notable ex- 



amples are: 



Barger et al 



yman 



2000 



break galaxies (L BGs; ISteidel et al.ll2003l ). sub -mm galaxies (SMGs; e.g. 



B 



Sawicki &: Thompson 



ain et al.ll2002h . Lya emitters (LAE gj e.g.. iFvnbo et al.l l2003bl 12005 



2006 ), distant red galaxies (DRG s; iDaddi et al.l 12004 ). damped Lya 



absor bers (DLAs; [Mgller et al.l 120021 ; I Wolfe et al.l 120051 ) , and high-redshift dropout galaxies 



(e.g., iBouwens et al.ll2010l ). Studies of these populations have led to major advances in our 
understanding of galaxy evolution, but our view remains far from complete, particularly at 
higher redshifts. 

The major stumbling blocks to the routine use of GRBs as tools for the investigation 
of cosmic history have been the relatively slow buildup of statistically useful samples, and 
the disconnect between the high-energy properties of the burst and parameters useful for 
measuring the host and its star formation rate (SFR). These latter parameters, such as the 
redshift, the host-galaxy bolometric luminosity, metallicity, and dust content can only be 
determined from longer wavelength observations in the UV, optical, infrared, millimeter, 
and radio regimes. 

Accurate X-ray localizations, as an intermediate step between the gamma-rays and the 
optical/infrared, have turned out to be critical to finding th e afterglows and th e host galaxies 
of GRBs. The first detections provided by Bepp oSAX ( jCosta et al.l Il997| ) heralded the 
afterglow era. Since 2005, Swift (IGehrels et al.ll2004f) : with its autonomous slewing capability 
and onboard X-Ray Telescope (XRT, iBurrows et al.ll2005l ) has been extremely successful in 
this regard, delivering high-accuracy (few arcseconds) positions for approximately a hundred 
bursts every year. It has become the GRB facility par excellence thanks to its innovative, 
intelligent design, open access policy, energetic team, and active community support. 

In spite of Swift' 's success, some of the problems that beset most other star-formation 
indicators also affect GRBs. Despite man y attempts to infer distances from the prompt 
high-energy properties (e.g.. lGuidorzill2005l ). UV/optical/NIR observations (and preferably 
spectroscopy) remain the only route to reliably find the redshift of a GRB or its host, and 
we must observe the hosts of GRBs directly across a range of wavelengths to connect SFR 
to GRB rate. Thus, to avoid optical biases in GRB selection requires a sample as com- 
plete as possible in redshift and host-galaxy properties. Swift's excellent X-ray localization 
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completeness and large sample size provide a first opportunity for this to be attempted. 



Several previo us studies of GRB hosts have found that the majority of GRB hos ts are 



faint, blue galaxies 



Fruchter et al. 



1999; 



Le Floc'h et al. 



low stellar masses rtCastro Ceron et all bood bold : ISavaglio et aJbopgh. T he GRBs gen 



20031: Ichristensen et al.lbo04h with 



erally occur within the UV-bright parts of their hosts (IBloom et al. 



2002: Fruchter et al. 



20061 ). consistent with their association with star formation. They are also typically fainter 



and have lower metallicities than the galaxies selected in Lyman-break su rveys in the same 



redshift range (IJakobsson et al.l 120051 ; ISavagliol 120061 ; iFynbo et al.l l2008bl ) . Accounting for 



lem ( 


Fvnbo et al. 12001. '. 


2009: De Pasaua 


— — „ — — , ~. 
le et al. 2003: Jal 


— — ~ — ^ — 

cobsson et al. 2004; Perlev et al. 


2009; 


Kriihler et al. 


2011 


Svensson et al. 


2012: 


Rossi et al. 


2012: 


Mel 


andri et al. 


2012' 


, has 


also been attempted with some success (e.g.. Jakobsson et al. 


2006b: 


Fvnbo et al 


200< 


)). 



Submillimeter (submm) and radio s urveys of GRB hos ts have shown that only a few 
could be ultraluminous infrared galaxies (IBerger et al.l 120031). indicating t hat obscured star 
formation is not dominant in GRB-selected galaxies (ITanvir et al.l l2004j ) . Moreover, the 
candidate s ubmm bright GRB hosts are bluer, hotter, a nd at lower redshift than typical of 



SMGs (e.g., IBerger et al.l 120031 : Michalowski et al.l 12008( 1 



Most of the studies mentioned above have been limited to fairly low-redshift samples 
(z < 1.5), which calls for an extens ion of such statistic al investigations to higher redshifts. 
Recently, more complete surveys (IGreiner et al.l l201ll ) have indicated the existence of a 
fractio n of GRB host gal axies being more chemically evolved, more dusty and of higher 
mass (IKruhler et al.l 1201 ll ). This highlights the importance of minimizing optical selection 
effects in usi ng GRBs as tracers of high-redshift star fo r mation and the cosmic s t ar-formation 



histo r y (e.g.. ITanvir et al. 



2012 



Basa et al. 



20121 ; 



2012 



Trent i et al. 



20121 ). 



Salvaterra et al. 2012: Robertson fc Ellis 2012; Elliott et al. 



We have conducted a survey of galaxies hosting GRBs, suitable for statistical studies 
and for observations with other existing or future facilities. In particular, we have imaged 
69 GRB host galaxies to deep limits in R and K s and obtained follow-up spectroscopy to 
determine redshifts and study Lya emission in the higher-redshift systems. The improved 
redshift completeness is also important for statistical studies of GRB properties themselves. 
The details of the survey design and target selection are described in Section 2. The observa- 
tions are described in Section 3, along with the catalog of GRB and host galaxy properties, 
including magnitudes and redshifts. Section 4 provides a summary of some of the immediate 
scientific results and Section 5 concludes. More detailed discussion and resu lts are presented 



i n fou r companion papers by D. M alesa ni et al. (2012 , in p reparation), IJakobsson et al. 



( 120121 ). iMilvang- Jensen et al.l (120121 ). and Kriihler et al.l ( 120121 ). Supplementary radio, Hub- 
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hie Space Telescope (HST), and X-ray observations are reported in iMichalowski et al.l ( 120121 ). 
S. Schulze et al. (2012, in preparation), and D. Watson et al. (2012, in preparation). 

All magnitudes given in this paper are in the Vega system. Throughout our work, we 
adopt a cosmology with Q m = 0.27, fl^ = 0.73, H = 71 km s _1 Mpc -1 . 



2. Target selection and sample definition 



A statistically meaningful and optically unbiased sample is essential in addressing issues 
such as the GRB redshift distribution, the importance of d ust obscuration, and th e determi- 
nation of the burst and host luminosity functions (see also ISalvaterra et al.ll2012l ). We have 
applied the following selection criteria to provide a large and homogeneous sample of targets 
and so as to ensure that selection effects are minimized and well understood. The criteria 



are akin to those introduced by lJakobsson et al.l (j2006bl ) and refined by lFynbo et al.l (120091 ). 



1. Our starting point is the sample of Swift-detected GRBs. Swift is the only current 
facility which can provide a large, homogeneous, well-localized sample of GRBs. We 
pick only events triggered onboard by the Burst Alert Telescope (BAT) and disregard 
those (very few) discovered during ground analysis. 

2. We limit ourselves to lo ns-duration GRBs, and discard all events wit h T 90 < 2 s 
( Kouveliotou et aDll993l ). based on the catalog of Sakamoto et al. ( 2011 ). 



The bursts must have a detected X-ray afterglow, the location of which is made avail- 
able less than 12 hr after the trigger. The requirement of the existence of an X-ray 
afterglow effectively makes our sample X-ray selected. 



The above criteria are essentially constraints on the nature of the GRBs. The require- 
ment on the timing in Criterion 3 additionally ensures that efficient afterglow searches could 
be performed. It also works as an effective Sun constraint, since those bursts with lines of 
sight close to the Sun upon discovery cannot be observed quickly. It should be noted that 
the fraction of long GRBs without a detected X-ray afterglow is very low, providing the XRT 
is repointed rapidly (see below). 

Looking for the hosts of all these GRBs would be very demanding. We therefore apply 
further cuts, with the purpose of reducing the sample size in an unbiased way. Imposing 
restrictions 4-6 does not bias the sample toward optically bright afterglows; instead each 
GRB in the sample has favorable observing conditions, i.e., useful follow-up observations are 
likely to be secured. 
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4. The Galactic foreground optical extinction as determined by lSchlegel et al.l (119981 ) must 
be Ay ^ 0.5 mag. This makes the detection level almost uniform (i.e., no bias against 
systems that happen to be faint just because of Galactic extinction). As a bonus, all 
fields have high Galactic latitude > 15°) ensuring that crowding by foreground 
stars is not an issue (Figured]). 

5. As another visibility constraint, we require that the distance on the sky between the 



GRB and the Sun at the time of explosion is more than 55° (similar to lJakobsson et al. 



2006bj). This is relevant since bursts closer to the Sun than this will offer very limited 
opportunity for ground follow-up. In practice, it is not a very demanding constraint 
because Criterion 3 effectively enforces a solar angle constraint of this order to satisfy 
the Swift Sun-avoidance condition for rapid slews for the narrow-field instruments 
(XRT and the Ultraviolet-Optical Telescope, UVOT) of 47°. In fact, during the survey 
period, only GRB070810A had an X-ray afterglow announced within 12 hr but had a 
discovery line of sight closer to the Sun. 

6. We require that there is no nearby star or bright contaminating object that would 
render host-galaxy identification and spectroscopy difficult. To make this quantitative, 
we note that, under the typical observing conditions of our i?-band imaging survey, 
bright stars affect a region of radius r ~ 1.8 + 0. 4 x exp[(20 — _R)/2.05] arcsec, where 



R is the USNO-B magnitude ( iMonet et al.l 120031 ) . This relation is valid in the range 



11 ^ R ^ 19- We thus require that no star with R < 19 is located on the CCD detected 
within this distance from the X-ray error circle edge. In two cases, GRBs 060604 and 
070721B, very bright stars (R < 10) at distances > V would nominally violate the 
above limit, but their effect was minimized by offsetting the pointing such that they 
were located outside of the detector. 



The aforementioned six criteria are applicable in principle to any survey of Swift GRB 
positions. Our survey was carried out with the VLT, and was restricted to bursts detected 
during a specific 2.5 yr period, so that for this sample we applied two additional criteria. 

7. Bursts triggered between 2005 March 1 and 2007 August 10. During this period Swift 
was fully operational and automatic slews were routinely enabled. There were 236 
GRBs (229 triggered) in this time@ 



2 We do not consider GRBs 060204C, 060602B, and 070610, which were very likely due t o Galactic sources 
(|Sakamoto et allbood: IPalmer et allbood iKasliwal et~al]|2008t ICastro-Tirado et alJliooj ). 
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8. In order to be well placed for observation with the VLT we limit the declination range 
to be -70° < 5 < +27° (J2000.0). 

These criteria can easily be extended by using other observatories and different points 
in time. 

After enforcing these criteria, we end up with a sample of 69 GRBs. The advantage 
of our selection is illustrated by the large success of afterglow follow-up for bursts in this 
subset. Based on Jochen Greiner's Web siteJU optical/NIR afterglows have been discovered 
for 75% (52/69) of the bursts in our sample, compared to 66% for the overall Swift sample. 
More importantly, spectroscopic redshifts are available for 55% of our targets (38 GRBs) 
compared to only 37% in the overall Swift sample. 

Our final criterion concerns the localization accuracy of the GRBs. 

9. The X-ray localization uncertainty must be better than or equal to 2 arcsec (90% 
error radius). This constraint is dictated by the necessity to perform meaningful host 
searches. Larger error boxes would lead to too many ambiguous detections. For bursts 
without an optical/NIR localization, w e ad opt the best constr ained position among 



the catalogs provided by iButlerl (l2007bl ) and lEvans et al.l (120091 ). In principle, a more 
accurate localization could be available from optical or radio data, but we only base our 
selection on the X-ray position, to retain a dependence solely on high-energy properties. 

We summarize the selection criteria in Table 1. Our final sample consists of 69 GRBs 
and is presented in Table 2. Figure [1] shows the distribution in the sky of all GRBs detected 
by Swift during the chosen time period with the TOUGH sample explicitly marked. 

Below we discuss the effect s of the various selecti ons. Figure [2] shows the 1 s BAT 15- 



150 keV peak photon flux (from ISakamoto et al.ll201ll ) versus the date. As can be seen, our 



targets uniformly sample the full BAT distribution and there is no obvious long-term change 
in the detection threshold as a function of time, which seems to be consistently at a level 
of about 0.4 7 cur 2 s~ l (the faintest target GRB 060218 is detected at 0.25 7 cm -2 s _1 ). 
This shows that we are starting from a stable, well-controlled set of events (Criterion 1), 
although we note that the effective trigger threshold does vary somewhat with many factors, 
such as the current background level, the position of the burst within the BAT field, and 
also the form of the particular GRB light curve (for example, the detection of some bursts 



3 http : / / www . mpe . mpg . de/~ j eg/ grbgen . html 
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is only made after integrating up more than the 1 s binning we consider here). Out of the 
69 TOUGH bursts, 56 (13) were rate (image) triggers. 

Our Criterion 3 is the detection of an X-ray afterglow. During the considered period, 
only 25 triggered GRBs were not detected (or observed) by XRT (11% of the overall Swift 
sample). Restricting to long-duration bursts (Criterion 2) with prompt slew (less than 
3600s), however, we find that only 5 events out of 180 (2.8%) were not detected by XRT. 
This shows that, in practice, XRT selection has a negligible impact on the size of our sample, 
although it does mean that rare GRBs with very faint X-ray afterglows will be systematically 
underrepresented. 

Figures [3] and H] illustrate the distributions of the times to first XRT observation (taken 
from the Swift GRB tabkjl; akin to but not identical to Criterion 3) and the Galactic V- 
band extinction ( ISchlegel et al.lll998l Criterion 4) for the TOUGH sample as well as the full 
distribution of GRBs. 

Criterion 9 introduces a potential bias in our requirement of an accurate XRT position, 
which might lead us to favor optically bright afterglows. The accuracy of the X-ray position 
is influenced by many factors, several of which are not directly relate d to the GRB properties. 
For example, the astrometrically corrected positions ( lButlerll2007bl ) require the detection of 
serendipitous X-ray sources in the field other than the GRB. This depends more on the 
total exposure time than on the afterglow brightness. Figure [5] shows that the accuracy 
of X-ray positions has been stable during the considered time span. In fact, many of the 
UVOT-enhanced positions are dominated by the systematic error term which is independent 
of the afterglow brightness. There is a possible dearth of very accurate positions after 2007 
March, which may be due to shallower exposures following a change of the Swift pointing 
strategy. In Figure El we investigate the impact of rejecting bursts without an accurate 
localization. We pick all GRBs which would pass sel ection criteria 1-8, regardless of the 
X-ray accuracy. We then plot their 90% error radii (jButlerl l2007al ; lEvans et al.l 120091 ) as 
a function of the cou nt rate at 12 hr, o btained by interpolating or extrapolating the light 
curves posted online§ ( Evans et al. 2009 ). Figure [6] shows that, indeed, faint afterglows have 
on the average larger error radii, but the effect is very small (a factor of ~ 3 across more 
than two orders of magnitude in the count rate). Thus, among 71 GRBs considered here, 
only 2 are discarded (less than 3%), showing that the impact on our sample is limited. 



4 http : / / swift . gsfc.nasa.gov/docs/ swif t/archive/grb_table/ 
s http : //www. swift . ac . uk/xrt .curves/ 



-9 - 



3. Survey and catalogs 



Our observational program consists of deep R and K s imaging and spectroscopy ob- 
tained with the F0RS1, F0RS2, ISAAC, and X-shooter instruments at the ESO VLT. 
All VLT observations were conducted at least 50 days after the GRB explosion in order 
to avoid possible contamination from the afterglow or an associated supernova. The pos- 
sible afterglow contamination in our images is quantified in D. Malesani et al. (2012, in 
preparation). Here, we briefly describe the overall results, as summarized in Table 3. De- 
tailed survev results are reported in four accompanying papers. Results on individual tar- 



gets or preliminary results are reported in '. 


rluiz-Velasco et al. (12007), Jakobsson et a . 


(2007. 


20091 2008 


2011a 


tJ3) 


Tanvir et al. 


(2008) 


, Thone et al. 


(20081. : 


20icf),|jaunsen et al. 


(2008), 


Xu et al. 


( 


2009|), 


Eh'asdottir et al. ( 


2009h. 


Ferrero et al. 


(2009), 


Fvnbo et al. 


(20091. : 


2008a), 


Malesani et al. (2 


!009, 


2010). Svensson et al. (2012). and 


Schulze et al. (2012 


). In Tables 2 



and 3, we also list 10 extra hosts that were observed for our program, but that do not obey 
all of the TOUGH sample criteria. 



3.1. Optical imaging 

For the purposes of identification of hosts and brightness estimation we have obtained 
deep images of the targets in the R band (typically rest-frame UV) to magnitude limits in 
the range R = 26-28 mag. The seeing FWHM of the i?-band images ranges from 0'.'50 to 
l'/26 with a median of 0'.'71. D. Malesani et al. (2012, in preparation) discuss our procedure 
for identifying a host galaxy, including a quantification of the chance probability of wrongly 
identifying an unrelated galaxy in the error circle as the host. We note that there is always 
a small risk of wrong host identification in both imaging and spectroscopy. While this could 
lead to erroneous conclusions when discussing individual systems, it will have little effect on 
statistical results inferred from a sufficiently large sample, like the one presented here, as 
long as the probability for false identification is low. 

We have identified host galaxies for 55 systems, corresponding to 80% of the sample. 
The coordinates of the identified host galaxies are reported in Table 3. In Figure [TJ we show 
a mosaic of GRB host galaxies that are only localized by the XRT (a mosaic of all identified 
host galaxies is presented in D. Malesani et al. (2012, in preparation). These 17 galaxies (25% 
of our sample) are typically missed from other surveys. The observed magnitudes are shown 
in Figure [HJ It is evident that host galaxies are generally sub-luminous, in the sense of being 
fainter than L* at z = 0. We note that the simple detection of a host ga laxy in the R band 



can be used to set an upper limit to the redshift of the GRB of z < 5.6 (IRuiz-Velasco et al. 



20071 : IJakobsson et al.ll2012f ) 
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3.2. Near-infrared imaging 



To obtain basic color information for the targets we also obtained iT s -band imaging 
(typically rest-frame visual) to a limiting magnitude of about 22. The seeing FWHM of the 
i^-band images ranges from / /34 to 0'.'83 with a median of 0?51. These data allow detection 
of extremely red objects (EROs, R — K s > 5) for all .R-band-detected targets. Although the 
majority of previously studied GRB host s have had blue op t ical-IR colors, a small number 
of ERO hosts have been identified (e.g . . iLeyan et al.l l2006at iBerger et al.l 120071 ; iHunt et al. 



2011: Svensson et al. 2012; Rossi et al. 



20121 ) . The resulti ng R — K x colors for o ur sample 



are shown in Figure El These are roughly consistent with iLe Floc'h et al.l (120031 ) . All host 
galaxies detected in K s were also detected in R. We have identified two EROs among our 
sample, the host galaxies of GRB050714B and GRB 070808 which both have R- K s ^5. 



3.3. Spectroscopy 



The identification of the hosts was used to determine their suitability for spectroscopic 
follow-up. For targets with no afterglow absorption redshift that are sufficiently bright for 
spectroscopic redshift determination (host galaxies generally are emission-line galaxies) we 
have obtained spectra with the purpose of making the redshift determinations for our sample 
as complete as possible. This is a crucial step in making any statistical inferences about the 
redshift distribution of GRBs and their host galaxies, in particular with reference to the star 
formation history of the universe. However, we note that the selection of galaxies at this 
stage based on their apparent .R-band magnitudes could introduce an optical bias into our 
sample, and that securing red s hifts for th e remaining (fa i nter) galaxies remains a priority. As 
reported by 1 Jakobsson et al.l ( 120121 ) and iKriihler et al.l (120121 ). we have obtained spectra of 
23 host galaxies in the TOUGH sample for which we have determined 15 new redshifts (and, 
unfortunately, found that 3 previously reported redshifts of the sample were erroneous). In 
total, we now have redshifts for 53 galaxies out of 69 in the sample. Moreover, even in cases 
where spectroscopy was unsuccessful in determining a redshift, spectra and imaging were 
used to set useful constraints, e.g., from an upper limit of the position of the Lyman limit 
or Lya breaks or from the absence of emission lines. 



4. Discussion 



We have verified that there are no obvious trends in the redshifts or magnitudes of the 
TOUGH host galaxies as a function of GRB epoch. 
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Below we summarize some of the results of the TOUGH survey, with special emphasis 
on the implications of having compiled an X-ray selected sample of GRB host galaxies, and 
obtaining photometry and spectroscopy of a sample of GRB host galaxies extending to faint 
and high-redshift systems. 



4.1. Redshift distribution 



In Figure 0, we show the redshift distribution prior to our survey and compare to the 
redshift distribution resulting from our survey. All redshifts are given in Table 3. The new 
redshift distribution is more smooth with a conspicuous peak at around redshift 2 with a 
median redshift of 2.14 ± 0.18. The median redshift of the entire TOUGH sample (i.e., 
including systems with no reds hift constraint) i s robu stly determined to be between 1.5 and 
2.5 (Table 4). As discussed by lJakobsson et al.l (120121 ) our redshifts and redshift constraints 
allow us to distinguish models for the GRB rate and luminosity function and formation 
biases. 

To illustrate the importance of sample selection criteria for the determination of the 
median redshift for GRBs we plot the redshift completeness as a function of BAT (15- 
150 keV) peak photon flux in Figure [TD1 The median 1 s peak flux for the TOUGH sample is 
1.4 7 cm -2 s -1 (Figure [TUT) . Subsamples with or without redshift give a consistent value for 
the median 1 s peak flux, as do subsamples with or without an optical/NIR afterglow. About 
75% of the GRBs in o ur sample have lower peak photon flux than the cutoff adopted by 



Salvaterra et al.l ( 120121 ) (> 2.6 7 cm 2 s 1 ). Our redshift completeness is in the range 70%- 



80% below that value and reaches almost 90% slightly above. The median redshift increases 
almost monotonically with lower peak flux, from 1.38 at a cutoff at 2.6 7 cm" 2 s -1 to 2.12 



for a cutoff at 0.4 7 cm s . Tying to the BATSE luminosity function, iSalvaterra et al. 
(120121 ) modeled their results and predicted a median redshift of z = 2.05 ±0.15 to that limit, 
consistent with our value. 

However, obtaining redshifts of fa i nter g alaxies may shift the median redshift somewhat 
higher. As indicated by lKriihler et al.l (120121 ) the faint host galaxies targeted with X-shooter 
have been found at significantly higher redshifts compared to brighter targets. We therefore 
suspect that the remaining very faint host galaxies may lead to even higher average redshifts. 
Figure [TTJ shows the median redshift as a function of the median i?-band magnitude of the 
hosts (the corresponding values are given in Table 4). The median redshift of the TOUGH 
sub-samples is lowest for the bright hosts for which the redshift is determined from emission 
lines using FORS, higher for fainter hosts where the redshift is determined using X-shooter, 
and higher still for bursts where the host remains undetected, ^-detected host galaxies are 
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at significantly lower redshifts than galaxies which are not detected in K s . The overall trend 
seems to be consistent with the expectation that fainter GRB host galaxies are on average 
at higher redshifts. 

TOUGH fo r the first time pic ks up high-redshift GRB ho sts in significant numbers. The 
compilations of I Chen et al.l ( 120091 ) and ISavaglio et al.l (120091 ) count eight optically detected 
GRB host galaxies at z > 2 of which two are Swift GRBs. Table 3 lists a total of 18 TOUGH- 
detected GRB host galaxies at z > 2, of which 17 are new, as well as one non-TOUGH z > 2 
GRB host galaxy. Such a large sample opens the field for more detailed studies of individual 
systems and of high-redshift GRB hosts as a whole (see also Section 14731 . 



4.2. Galaxy colors and magnitudes 



Of special interest is the study of dark bursts which are believed to be caused mostly by 



dust obscuration (e.g., 


Dioreovski et a 


2001; 


2008; 


Perlev et al. 


2009; 


Kriihler et al. 


2011; 



Rol et aJj 



Svensson et al 



2007kljaunsen et alihopj: 



201 



2: 



Rossi et al 



Tanvir et al 



20j3). Hence. 



it is expected that the properties of the host galaxies of XRT-only GRBs may be different 
from those hosting GRBs with detected optical/NIR afterglows. 



Figure [T21 highlights the difference between systems with and without an optical/NIR 
afterglow. There is a significant dearth of XRT-only galaxies at very faint magnitudes 
in the R band. This is also reflected in the median i?-band magnitude of 25.8 ± 0.3 for 
galaxies hosting an optically/NIR-detected afterglow, as opposed to 24.4 ± 0.4 for those 
with undetected ones (Table 4). XRT-only localised host galaxies are also relatively brighter 
in K s (median magnitude of K s = 21.1 ± 0.3 opposed to K s > 22.5 for galaxies hosting 
GRBs with optical/NIR afterglows). A clue to the origin of these remarkable differences 
comes from the lower panel of Figure [12] which demonstrates that XRT-only hosts have 
equi yalent X-ray absorp tion column densities at z = extending to much higher values (see 
also lKruhler et al.ll2012l ). These differences clearly demonstrate that statistical studies based 
on optically selected GRBs will be strongly biased. 

We note that our procedure of identifying the host as the brightest galaxy inside the XRT 
error circle naturally will lead to a bias toward brighter galaxies if any of the identifications 
are incorrect. To quantify this bias we tested our host identification procedure in the XRT 
error circles of GRBs with an optical/NIR afterglow. We found that in 6 out of 50 cases 
(12%) we would misidentify the host (this would indicate that about 1-2 of the galaxies 
in Figure [7] may have been misidentified). This procedure led to a biased median i?-band 
magnitude of 25.5 ± 0.2 for the XRT error circles of GRB with optical/NIR afterglows but 
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XRT-only identified host galaxies. Hence, the difference between host galaxies with and 
without optical/NIR afterglows remains significant. 

Regarding systems with both R and K s detections, we find that the median magnitude 
of hosts of GRBs with or without optical/NIR afterglows is consistently R ~ 24.3 ± 0.3. 
However, the median colors of XRT-only GRB host galaxies detected in the K s band are 
significantly redder, R — K s = 3.7 ± 0.3, than those for which an optical/NIR afterglow was 
detected, R — K s = 2.8 ± 0.1. While this difference in color could partly be a redshift effect, 
their larger X-ray absorption column densities suggest that GRBs without an optical/NIR 
afterglow suffer stronger attenuation due to dust and that dust extinction causes optical/NIR 
GRB afterglows to drop out of detection. This is in line with other studies which find that 
host galaxies of dar k GRBs are redder and brighter (more ma ssive) than those of GRBs with 
optical afterglows (jKruhler et ajj|201l| ; ISvensson et al.ll2012l ). A more detailed discussion of 
the R — K s distribution as a function of redshift and a comparison to other high-redshift 
galaxies are presented in D. Malesani (2012, in preparation). 



4.3. Star formation in GRB host galaxies 



Motivated by the fact that, before our sur vey, five out of a possible five hosts in a 



pre- S wift GRB sample were found to be LAEs (iFynbo et al.l 120021 . 12003al ; iJakobsson et al. 
20051 ). we have also performed spectroscopic observations of systems with known redshift in 
the redshift ra nge 1.8-4.5 where the Lyq e mission line can be best observed. Remarkably, 
as reported by iMilvang- Jensen et al.l (120121 ). only 7 out of 20 galaxies targeted have secure 
detection of Lya emission while upper limits on the Lya equivalent width are derived for 
another 7 galaxies from detecting the continuum in the spectra (Figure [S]). The origin of 
the lower (but more representative) success rate is likely a combination of lower intrinsic 
luminosities of the galaxies (and hence lower average Lya luminosities) and lower equivalent 
widths, at least partly due to higher intrinsic extinction in the host galaxy. The Lya emission 
is still stronger than in flux-limited samples of LBGs. Moreover, by comparison to the 
afterglow redshifts, we find that the velocity centroid of the Lya line is redshifted by 200- 
600 k m s" 1 with respect to the systemic velocity, indicative of outflows in the host galaxies 
(e.g., lAdelberger et al.ll2003l ; IVerhamme et al.ll2006l ; lLaursen et al.l 120091 ) . 



In order to assess the SFR s of GRB hosts, we have performed radio observations of 
TOUGH GRB hosts at z < 1 (IMichalowski et al.ll2012l ). This redshift limit was chosen 
to obtain meaningful limits on SFRs. We did not detect any TOUGH GRB hosts, which 
indicates that their average SFR is below ~ 15 M yr" 1 . We also found that at least 65% of 
GRB hosts at z < 1 have SFR < 100 M yr~ x and that at least 92% of them have A v < 3 
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mag. The obtained limits allowed us to conclude that the distribution of SFRs and dust 
attenuation of GRB hosts at z < 1 is consistent with that of other star-forming galaxies at 
similar redshifts. 



Finally, in S. Schulze et al. (2012, in preparation) the magnitudes are used to determine 
the luminosity function of GRB-selected galaxies. In particular, determining the slope of 
the faint end power-law part of the luminosity function is important in determi ning whether 



such galaxies constitute a major fraction of the star formation in the universe 


(Fvnbo et al. 


2001, 


2002; 


Jakobsson et al. 


2005. 


2012; 


Robertson & Ellis 


2012; 


Tanvir et al. 


2012). 



5. Summary 

Long-duration GRBs provide a complementary view of the nature of "typical" star- 
forming galaxies over most of cosmic history. However, to date, most samples of GRB 
hosts have been heterogeneous and had a strong bias toward GRBs with detected optical 
afterglows. In this paper we have described a large, homogeneous sample of hosts targeted 
with the VLT, using X-ray positions where no optical/NIR afterglow was found. We report 
the main catalog in Table 3. We detect the majority of the hosts, namely 80%, in the i?-band 
down to a detection limit of ~ 27.5, including 17 systems at z > 2. Within the TOUGH 
sample, we detect 54 hosts, of which 32 were not previously detected. GRB host galaxies 
cover a wide range of luminosities from < 0.1L* to L*. From the i^ s -band imaging we confirm 
the result of more biased studies that GRB host galaxies tend to be blue. We also find that 
galaxies with no optical/NIR afterglow are significantly brighter and redder than those with 
an optical/NIR afterglow, highlighting the importance of TOUGH in being X-ray selected. 
With our newly reported redshifts, the redshift distribution consisting of 77% of the GRBs, 
now appears smooth with a median redshift above 2. 

Our sample will form the basis for statistical studies of GRB host galaxies, such as their 
luminosity function, colors and SFRs, redshift distribution, Lya properties, radio emission, 
and X-ray absorption properties. 

Future observational work will focus on completing the sample, in particular in redshift, 
and to complement it at other wavelengths, in particular, radio, sub- mm (ALMA, Herschel), 
and NIR/MIR (Spitzer, HST, JWST). It is also useful for selecting subsamples for further 
study, such as SED studies or emission-line diagnostics, e.g., metallicity determination. Fi- 
nally, it is straightforward to enlarge in a well-defined way, in time or declination, as we are 
doing with VLT/X-shooter and Gemini. 
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Fig. 1. — All-sky map (Mollweid projection) of the 236 Swift GRBs which occurred between 
2005 March 1 and 2007 August 10 (Criterion 7). Filled circles: GRBs with Swift /XRT 
detection. Filled circles, encircled in red: GRBs obeying all our sample selection Criteria 
1-9. Empty circles: GRBs with no Swift /XRT detection (excluded by Criterion 3). Blue 
squares: GRBs classified as short (excluded by Criterion 2). Crosses: nontriggered GRBs 
(excluded by Criterion 1). The declination cuts are shown with dashed lines (—70° and +27°; 
Criterion 8) and the region with Galactic latitude |fe| < 15° is shown with the dotted curves, 
which roughly corresponds to the sample selection criterion Ay ^ 0.5 mag (Criterion 4). 
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Fig. 2. — 1-s Swift /BAT peak photon flux (15-150 keV) for all Swift GRBs observed during 
the time span of the TOUGH survey. Empty circles: long bursts without an X-ray afterglow. 
Filled circles: long GRBs with Swift /XRT detection. Filled, encircled circles: GRBs obeying 
all our sample selection Criteria 1-9. Squares: GRBs classified as short (Tg < 2 s). Crosses: 
nontriggered GRBs. There seems to be a constant and uniform detection lev el as a function of 



time (the very faint system at 0.04 7 cm s is the unusual GRB 060123; ICummings et al. 



20061 ). The vertical dashed lines bracket the s elected time span f or the sample. The horizontal 



dotted line is the peak flux limit adopted by ISalvaterra et al.l ( 120121 ) . 
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Fig. 3. — Distribution of times post-trigger to the first XRT observation for all Swift GRBs 
happening between 2005 March 1 and 2007 August 10. For 12 bursts, the XRT was not 
repointed. The fraction constituting the TOUGH sample is indicated in red. The vertical 
dashed line marks the sample limit t ^ 12 hr, which is fulfilled by 92% of the repointed 
bursts (207/224). Three GRBs (050603, 060223A, and 070810A) have £ X rt < 12 hr but the 
notice came later than 12 hr. These GRBs are indicated by arrows in the histogram. 
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Galactic K-band extinction (mag) 



Fig. 4. — Distribution of Galactic foreground extinction for all Swift GRBs detected between 
2005 March 1 and 2007 August 10 (13 bursts have 3 mag < Ay < 17 mag and are not shown 
in the figure). The fraction constituting the TOUGH sample is indicated in red. The 
vertical dashed line marks the sample limit Ay ^ 0.5 mag. Seventy-one percent of all bursts 
(167/236) have Ay sC 0.5 mag. 
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Fig. 5. — Swift/XRT error radius (90% encircled likelihood) as a function of date for GRBs 
detected during the time span of the TOUGH survey. Filled and empty points indicate 
bursts with and without a detected optica l/NIR afterglow (OA), respectively. Green squares 
indicate positi ons from the Butler d2007b) catalog, while blue diamonds are from the UVOT- 



enhanced list (IGoad et al. 



2007; 



Evans et al.ll2009l ). Our sample targets are circled in red. 



During the period 2005 March 1 through 2007 August 10 (marked by vertical dashed lines), 
89% of the GRBs have an XRT error circle radius ^ 2" (horizontal dotted line). The two 
dotted circles indicate the two bursts that are excluded only by Criterion 9 (r ^ 2"; see 
Figure |6]). The average error radius is l'!7. 
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Fig. 6. — X-ray error radius vs. the Swift /XRT count rate at 12 hr (0.3-10 keV). Only 
bursts which "survive" Criteria 1-8 are shown. Filled and empty points indicate bursts with 
and without a detected optical/NIR afterglow (OA), respectively. The horizontal dashed 
line indicates the TOUGH constraint set by Criterion 9. Notice the flat distribution, which 
implies that we only minimally bias our sample. Indeed, only two GRBs are excluded. The 
dotted line indicat es the systematic error contributing to the UVOT-enhanced uncertainties 
ffavans et alJboQflh . 
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Fig. 7. — Mosaic of the fields of 17 XRT-localized GRBs with no optical/NIR afterglow. 
Each panel, showing the i?-band image, is 15" x 1 5" in size (north is up a nd east is left) . 
Blue and green circles identify the XRT positions by Evans et al.l (120091 ) and iButlerl ( l2007bl ). 
respectively. 
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Fig. 8. — TOUGH sample properties as a function of redshift. Optically localized systems 
are shown as blue filled circles, XRT localized systems are shown as open red circles. Upper 
limits are shown as arrows. Hosts without a reported redshift are plotted on the left side of 
the diagram, (a) i?-band host magnitudes (Vega system). The dashed curve corresponds to 
the observed i?-band magnit ude for an Mr = —21 ( Vega) object assuming F v oc z/ -0 ' 5 and 
IGM opacity as described by iDalFAglio et al.l ( 120081 ). (b) .fr s -band host magnitudes (Vega 
system), (c) R — K s colors. The ERO limit of R — K s = 5 is shown as a dotted line. Tw o 
EROs are detected, (d) Lya detections and upper limit s from IMilvang- Jensen et al.l ( 120121 ) . 
The upper limit reported for GRB 05 0904 at z = 6.30 (ITotani et al.l 120061 ) was adjusted to 
be comparable with those reported in IMilvang- Jensen et al. d2012h . The Lya fluxes are in 
units of 10~ 18 erg cm -2 s _1 . The vertical dotted lines mark the redshift interval where Lya 
was searched in our data. 
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Fig. 9. — Redshift distribution of the TOUGH sample. In both panels the dashed vertical 
line indicates the median redshift. (a) Without considering the TOUGH observations, 38 
redshifts were considered secure with a median redshift of z — 2.35. (b) Our TOUGH spec- 
troscopic observations add 15 new redshifts and dem onstrate that three redshifts reported 
in the literature are erroneous ( iJakobsson et al.ll2012[ ). Here the median redshift is z = 2.14. 
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Fig. 10.— Redshift and BAT 15-150 keV p eak photon flux . The vertical dashed line 
indicates the peak photon flux limit adopted by ISalvaterra et all ( 120121 ) . (a) The cumulative 
distribution of peak photon fluxes of the TOUGH bursts, (b) The redshift completeness 
fraction of the TOUGH sample bursts as a function of the limiting peak photon flux. The 
redshift completeness is constant around 75% below ~ 3 7 cm~ 2 s _1 . In order to significantly 
increase the completeness fraction, one would have to discard at least > 75% of all bursts, (c) 
The median redshift of TOUGH GRBs with known redshift brighter than the corresponding 
peak flux. Fainter GRBs are on average at higher redshift. 
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Fig. 11. — Median redshift (with a range of uncertainty where the median is not uniquely 
determined) against median i?-band host magnitude for various subsets of our sample. In 
the case of those without an i?-band detection, the set has been placed at R ~ 27. The 
figure broadly illustrates that the fainter hosts are typically at higher redshift, characteristic 
of our magnitude limit. The subset of hosts without redshift has a median R w 25.9 which 
suggests a median redshift above 2 if it follows the same trend. 
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Fig. 12. — Distributions of host galaxy R- (top) and f^-band (middle) magnitudes as well as 
X-ray absorptions in excess of Galactic, determined at z = (bottom) for the TOUGH sam- 
ple (open histogram). Subsets of galaxies selected with a localization from an optical/NIR 
afterglow (OA) and a localization from Swift /XRT only are shown as hatched and solid 
histograms, respectively. The corresponding median values (or upper limit to the median 
as indicated by arrows) are shown as dashed or dotted lines. The histograms were created 
from the sum of the normalized probability distributions for each object derived from their 
uncertainties or upper limits. For the galaxy magnitudes the upper limit probabilities were 
evenly distributed in magnitude space from the upper limit to a lower cutoff value: R = 30 
and K s = 26. While the height of the histogram at low values depends on this choice, the 
histogram area does not. 



Table 1. Sample Selection Criteria 



Criterion No. Short Description Type of Constraint 



1 Onboard Swift /BAT trigger GRB properties 

2 Long GRBs: duration T90 > 2 s GRB properties 

3 X-ray afterglow position available within 12 hr GRB properties, visibility, and slight biasing 

4 Galactic Ay ^ 0.5 mag Visibility 

5 Sun distance > 55 deg Visibility 

6 Not near bright stars Visibility 

7 2005 Mar 1 - 2007 Aug 10 VLT sample 

8 -70° < 8 < +27° (J2000.0) VLT sample 

9 Localization accuracy better than 2" Slight biasing 
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1.4 


UVW1 


9 


0.04 


62 


10.0 


1.91 


060714 


219101 


15:11:26.46 


-06:33:58.8 


1.4 


j3 


37 


0.26 


99 


116.0 


1.27 


060719 


220020 


01:13:43.70 


-48:22:50.6 


1.4 


if 


38 


0.03 


77 


66.9 


2.21 


060729 


221755 


06:21:31.79 


-62:22:12.5 


1.4 


UVW2 


39 


0.18 


124 


113.0 


1.14 


060805A 


222683 


14:43:43.47 


+12:35:11.2 


1.6 


X 


40 


0.08 


93 


4.9 


0.31 


060814 


224552 


14:45:21.36 


+20:35:09.2 


1.4 


if 


41 


0.13 


72 


144.9 


7.26 


060908 


228581 


02:07:18.42 


+00:20:32.2 


1.4 


UVW1 


9 


0.10 


72 


18.8 


3.09 


060912A 


229185 


00:21:08.13 


+20:58:18.5 


1.4 


UWM2 


9 


0.17 


109 


5.0 


8.27 


060919 


230115 


18:27:41.80 


-51:00:52.1 


1.7 


X 


42 


0.24 


87 


9.0 


2.31 
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Table 2 — Continued 



GRB 


Trigger No. 


R.A. 


Decl. 


Error 


Band 


Reference 


A v 


*XRT 




Peak Flux 






(J2000.0) 


(J2000.0) 


(") 






(mag) 


(s) 


W 


(7 cm -2 s~ v ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


060923A 


230662 


16:58:28.14 


+12:21:37.9 


1.5 


K 


43 


0.20 


81 


51.5 


1.38 


060923B 


230702 


15:52:46.70 


-30:54:13.7 


1.8 


X 


■14 


0.49 


111 


8.9 


1.43 


060923C 


230711 


23:04:28.28 


+03:55:28.1 


3.5 


J 


45 


0.21 


203 


67.4 


0.89 


060927 


231362 


21:58:12.05 


+05:21:49.0 


1.6 


1 


46 


0.21 


65 


22.4 


2.68 


061004 


232339 


06:31:10.81 


—45:54:24.2 


1.4 


X 


47 


0.16 


60 


6.3 


2.48 


061007 


232683 


03:05:19.60 


-50:30:02.4 


1.4 


UVW2 


48 


0.07 


80 


75.7 


14.70 


061021 


234905 


09:40:36.12 


-21:57:05.2 


1.4 


UVW2 


9 


0.19 


78 


43.8 


6.12 


061110A 


238108 


22:25:09.89 


-02:15:30.4 


1.6 


R 


49 


0.30 


69 


44.5 


0.48 


061110B 


238174 


21:35:40.46 


+06:52:32.9 


1.7 


R 


50 


0.13 


3195 


132.8 


0.47 


061121 


239899 


09:48:54.59 


-13:11:42.1 


1.4 


UVW2 


51 


0.15 


55 


81.2 


21.30 


070103 


254532 


23:30:13.80 


+26:52:34.4 


1.5 


X 


52 


0.22 


69 


18.4 


1.05 


070110 


255445 


00:03:39.38 


—52:58:28.1 


1.4 


u 


53 


0.05 


93 


79.7 


0.62 


070129 


258408 


02:28:00.98 


+11:41:03.4 


1.4 


R 


54 


0.46 


134 


459.1 


0.56 


070224 


261880 


11:56:06.57 


-13:19:48.8 


1.6 


R 


55 


0.19 


143 


48.0 


0.34 


070306 


263361 


09:52:23.29 


+10:28:55.5 


1.4 


H 


56 


0.09 


153 


209.2 


4.04 


070318 


271019 


03:13:56.77 


—42:56:47.3 


1.4 


UVW2 


57 


0.06 


64 


131.5 


1.64 


070328 


272773 


04:20:27.68 


—34:04:00.7 


1.4 


X 


58 


0.12 


88 


72.1 


4.22 


070330 


273180 


17:58:10.20 


—63:47:34.4 


1.5 


V 


59 


0.21 


68 


6.6 


0.88 


070419B 


276212 


21:02:49.77 


-31:15:48.7 


1.5 


V 


60 


0.30 


81 


238.1 


1.38 


070506 


278693 


23:08:52.31 


+10:43:20.8 


1.8 


B 


61 


0.13 


127 


1.3 


0.95 


070611 


282003 


00:07:58.12 


-29:45:20.4 


1.8 


UVW1 


62 


0.04 


3300 


13.2 


0.82 


070621 


282808 


21:35:10.08 


-24:49:03.1 


1.4 


X 


63 


0.16 


111 


33.3 


2.50 


070721B 


285654 


02:12:32.93 


-02:11:39.7 


1.5 


V 


61 


0.10 


92 


336.9 


1.57 


070802 


286809 


02:27:35.78 


-55:31:39.5 


1.5 


B 


65 


0.09 


138 


15.8 


0.40 


070808 


287260 


00:27:03.36 


+01:10:34.4 


1.5 


X 


66 


0.08 


114 


58.4 


1.94 


050412 


114485 


12:04:25.28 


-01:12:00.6 


3.7 


X 


67 


0.07 


99 


26.5 


0.49 


050603 


131560 


02:39:56.89 


-25:10:54.9 


1.5 


V 


68 


0.09 


39024 


22.0 


15.30 


060117 


177666 


21:51:36.23 


-59:58:39.3 


1.5 


R 


69 


0.12 




16.9 


48.80 


060223A 


192059 


03:40:49.56 


-17:07:49.7 


1.5 


V 


9 


0.39 


86 


11.3 


1.35 


060313 


201487 


04:26:28.46 


-10:50:41.2 


1.1 


UVW2 


70 


0.23 


79 


0.7 


10.80 


060505 


208654 


22:07:03.32 


-27:48:53.0 


2.0 


B 


71 


0.07 


51744 


4.0 


2.65 


060602A 


213180 


09:58:16.28 


+00:18:14.4 


5.2 


R 


72 


0.08 


155520 


74.7 


0.56 


070721A 


285653 


00:12:39.13 


-28:33:00.9 


1.6 


X 


73 


0.05 


86 


3.4 


0.68 


070810A 


287364 


12:39:51.24 


+10:45:03.2 


1.4 


V 


74 


0.07 


88 


9.0 


1.89 


080207 


302728 


13:50:02.93 


+07:30:07.9 


1.4 


X 


75 


0.08 


121 


340.0 


1.00 



Note. — Colu mn 1: GRB name; Column 2: Swift trig ger number; Column s 3-5: X-ray position and error {Goad et al.ll2007l: 
lEvans et aljbooih : the position of GRB 060117 is taken from ljelmek et all ll2006T ); Column 6: bluest optical/NIR afterglow detection, 
excluding white and unfiltered measurem ents; "X" indicates n o optical/NIR detection; Column 7: reference for the afterglow detection; 
Column 8: Galactic V-ba nd absorption llSchlegel et al,lll998h ; Column 9: time to start of the XRT observatio n; Column 10: observe d 
gamma-ray duration Tgo llSakamoto et all201lT h" Column 11: peak photon flux over 1 s timescale (15-150 keV; ISakamoto et al.l201 ill . 
The horizontal line separates the o bjec ts in the TOUGH sa mple from 10 extra sy ste ms that were observe d wi thin our program. 
Ref erences. (lllBereer et all l (2005h; (2llStm et all j2005h; OllKamble et all d2009f); f4llSchadv et ail | |2006|1; (5l | Holland et all (120071) ; 

(11) 
(16) 



(61 lAfonso et al.N201lh; (7) Blustin et all ll2006ll ; (8) iLevan et aUfeOOdl ; (91 lOates et all |2009l); (101 IPandev et all ( 20Qi 
IPerri et all J2007h ; (121 Ide Pasquale et aUll2007l l; f 131 iBand et all d2005h : (14) iKennea et all ll2005ll ; (15l lSchadv et all I20T 
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iBlustin et al] [|2005al); (17)ISchadv et al.1 j2007bh; (18)lTagliaferri et al] d2005h ; (19 ) lde Pasquale et a.1.1 j2005h ; (20) lBloom 



200E ); (21)l Norris et all hood: (22) 



2005); (26) 



_ . __ . Moretti et 
^^^^ ^awa*^T*a L| l|2006l) 



27 



20051); (23) iMorris et a.1.1 (1200511; (24) iBlustin et al.l i2005bl); (25 



Alatalo 



Campana et al.l 
hold): 



hOOd): (28) 



(32) iBlustin 
ISchadv fc Morettil ([2006 ) ; 



Angclini ct al 
Page! ||200&!); 



hood: 



(29] 
(200 



Bcardmorc ct al. 
D'Avanzo ct al. I hood: 



(33) iFerrero et alJ (I2009D; (34)~lziaeepour et al.l 
Asfandvarov et al.l hood): (3 8) iMalesani et al.l J2006I 



30) iDai et al.l I 20071); (31) iThone et al 
l200Sft; (35) iMangano et al.l J2007f); (36] 

(39) | Grupe et al. (12007); (40) jZiaeepour etafl [|2006bl); (41) jLevan et al.l d2006bh; (42) iGuidorzi et al.l hood: (43) iTanvir et al 
hOOS): (44)IStamatikos et al.l ||2006|); (45) ICovino et al.M2006h; (46) iRuiz-Wasco et al.1 COOjl); (47) IZiaeepour et al.N2006al); (48 ) 



Sakamoto et al 



20071); 



(57) IChester et a.1.1 



ISchadv et al.1 J2007ah; (49) IChen et all I200I ); (50) iMelandri et all feood); (51) IPage et al.M2007l); ( 5 2) ISal 
(53) iTroia et al.l d2007h; (54) IMalesani et al. ll2007l)r^5) iThoene et al.l d2007al); (56) Ijaunsen et all ll2008h. 

hood): (58)lMarkwardt et all d2007l); (59)lKuin et a,U j2007i) ; ( 60)lde Ugarte Postigo et al.[ ||2007| ); (61 ) bandsman fc Paganil |2007l) ; 
(62) lLandsman et al.l J2007h ; (63) ISbarufatti et al.l ll2007^(64) lde Pasquale fc Ziaeepoud I2OO7I); (6 5) lElfasdottir et al.l (120091) ; (66) 
Cummings et all ||2007|) ; (67)lMineo et al.N2007h; (68)lGrupe et al.l I 20061) ; ( 69) I Jelmek et al.l \200d : (7 0)lRoming et al l \200d : (71) 
Xu et alJ j2009h ; (72) Ijensen et al.l hood : (73) Iziaeepour et al.l i2007l ); (74) IChester et al.l J2007h ; (75) iRacusin et al.l i2008l ). 
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Table 3. Catalog of Survey Results 



GRB R.A. Decl. R 

(J2000.0) (J2000.0) (mag) 

(1) (2) (3) (4) 



K s Redshift Notes Reference 

(mag) 

(5) (6) (7) (8) 



050315 


20 


:25 


:54. 


.189 


-42: 


:36 


:02 


16 


24.43+0. 


15 


21.89+0.36 


1.9500 


A 


1 


050318 


















> 26.79 




> 22.31 


1.4436 


A 


1 


050401 


16 


:31 


:28. 


794 


+02. 


11 


:14 


.33 


26.15+0 


31 


> 21.73 


2.8983 


A* 


2,3 


050406 


02 


:17 


:52. 


224 


-50: 


ill: 


:14 


.97 


26.61+0. 


,34 


> 21.25 


2.70 


P 


4 


050416A 


12 


:33 


:54. 


602 


+21: 


:03: 


26 


.58 


23.13+0 


.03 


20.36+0.20 


0.6528 


E 


1,5 


050502B 


















> 26.79 




> 21.81 


5.2 


P 


6 


050525A 


18 


:32 


:32. 


590 


+26: 


:20 


:22 


.60 


25.95+0. 


,34 


> 21.23 


0.606 


AE 


7 


050714B 


11: 


:18 


:47. 


717 


-15: 


:32: 


:49 


01 


25.45+0 


.20 


20.39+0.17 








050726 


















> 26.14 




> 21.28 








050730 


















> 27.20 




> 21.36 


3.96855 


A* 


8,9,3,10 


050801 


















> 26.78 




> 21.70 


1.38 


P 


11,12 


050803 


23 


:22 


:37. 


850 


+05: 


:47: 


:08 


90 


26.18+0 


33 


> 21.50 








050819 


23 


:55 


:01. 


598 


+24: 


51 


:38 


16 


24.10+0. 


.09 




2.5043 


E 


13 


050820A 


22 


:29 


:38. 


114 


+ 19: 


:33: 


:36 


61 


24.83+0. 


.18 


> 21.18 


2.61469 


A* 


14,15 


050822 


03 


:24 


:27. 


214 


-46: 


01 


:59 


43 


24.19+0 


.08 


20.78+0.20 


1.434 


E 


16 


050824 


00 


:48 


:56. 


151 


+22: 


:36 


:32 


.80 


24.23+0. 


.16 


> 21.41 


0.8278 


AE 


17,3 


050904 
























6.295 


A 


18 


050908 


01 


:21 


:50. 


730 


-12 


:57: 


17 


30 


27.76+0 


45 


> 21.91 


3.3467 


A* 


3,19 


050915A 


05 


:26 


:44 


843 


-28: 


:00 


:59 


.62 


24.56+0. 


.16 


20.69+0.24 


2.5273 


E 


13 


050922B 


00 


:23 


:13. 


345 


-05: 


.36 


:18 


31 


26.18+0. 


.30 


> 22.17 








050922C 


















> 26.34 




> 21.47 


2.1992 


A* 


20,21 


051001 


23 


:23 


:48, 


800 


-31: 


:31: 


23 


45 


24.36+0 


13 


21.13+0.35 


2.4296 


E 


13 


051006 


07 


:23 


:14. 


109 


+09: 


:30 


:19 


.96 


22.99+0. 


.07 




1.059 


E 


16 


051016B 


08. 


:48 


:27. 


844 


+ 13: 


:39: 


20 


.07 


22.64+0 


.03 


20.04+0.20 


0.9364 


E 


22 


051117B 


05 


:40 


:43, 


287 


-19: 


16: 


26 


.31 


21.06+0 


02 


17.75+0.05 


0.481 


E 


16 


060115 


03 


:36 


:08. 


355 


+ 17: 


20 


:42 


.65 


27.14+0. 


.53 


> 21.27 


3.5328 


A* 


3 


060218 


03 


:21 


:39 


683 


+ 16: 


:52: 


:02 


.09 


19.98+0. 


.01 


17.94+0.09 


0.03351 


E 


23,24,25,26,27 


060306 


02 


:44 


:22. 


930 


-02: 


:08: 


:54. 


,55 


24.09+0 


.08 


19.95+0.17 








060522 


















> 26.64 




> 23.15 


5.110 


A 


28 


060526 


















> 27.02 




> 22.11 


3.2213 


A 


29,20,30 


060604 


22 


:28 


:55. 


.029 


-10: 


:54: 


:55. 


.72 


25.52+0 


.18 


> 21.37 


2.1357 


AE 


3,13 


060605 


















> 26.54 




> 21.12 


3.773 


A 


31 


060607A 


















> 27.92 




> 21.86 


3.0749 


A 


15 


060614 


21 


:23 


:32. 


103 


-53: 


:01 


:36 


.20 


22.53+0 


06 


20.36+0.20 


0.125 


E 


32,33 


060707 


23 


:48 


:19. 


061 


-17: 


:54: 


:17 


72 


24.86+0 


06 


> 19.71 


3.4240 


A 


20,3 


060708 


00 


:31 


:13. 


748 


-33: 


:45: 


:32 


.73 


26.76+0. 


,28 


> 21.04 


1.92 


P 


11 


060714 


15 


:11 


:26. 


436 


-06 


33 


:58 


.24 


26.45+0 


,28 


> 21.40 


2.7108 


AE* 


20,3 


060719 


01 


:13 


:43, 


706 


-48: 


:22: 


:51 


.31 


24.62+0 


.12 


21.25+0.27 


1.5320 


E 


13 


060729 


06 


:21 


:31. 


794 


-62: 


.22: 


:12 


53 


23.71+0. 


.04 


21.20+0.31 


0.5428 


A 


3 


060805A 


14 


:43 


:43. 


381 


+ 12 


35 


:10 


.28 


23.48+0 


.03 


20.69+0.15 








060805A 


14 


:43 


:43. 


488 


+ 12: 


:35: 


:12 


56 


25.11+0 


,14 


21.66+0.38 








060814 


14 


:45 


:21. 


347 


+20: 


:35: 


:10 


,72 


22.85+0. 


.11 


20.06+0.10 


1.9229 


E 


16,13,34 


060908 


02 


:07 


:18. 


406 


+00: 


:20 


:31 


.55 


25.53+0. 


.18 


22.53+0.43 


1.8836 


AE 


35,3 


060912A 


00 


:21 


:08. 


127 


+20 


:58: 


:17 


66 


22.68+0 


04 


19.86+0.09 


0.937 


E 


36 
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Table 3 — Continued 



GRB 


R.A. 


Decl. 


R 


K s 


Redshift 


Notes 


Reference 




(J2000.0) 


(J2000.0) 


(mag) 


(mag) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


060919 


18:27:41.790 


-51:00:50.94 


25.79+0.26 


21.41+0.38 








060923A 


16:58:28.175 


+12:21:38.96 


26.07+0.24 


> 21.25 








060923B 


15:52:46.566 


-30:54:14.57 


24.34+0.16 


20.39+0.19 








060923B 


15:52:46.583 


-30:54:12.70 


24.40+0.16 


20.59+0.20 








060923C 


23:04:28.266 


+03:55:29.29 


25.45+0.18 


21.73+0.34 








060927 








> 21.76 


5.4636 


A 


37,3 


061004 






> 25.76 


> 21.24 








061007 


03:05:19.587 


-50:30:02.43 


24.40+0.17 


21.06+0.37 


1.2622 


AE 


38,3 


061021 


09:40:36.124 


-21:57:05.05 


24.44+0.06 


> 21.14 


0.3463 


AE 


16,3 


061110A 


22:25:09.843 


-02:15:31.12 


25.65+0.26 


> 21.41 


0.7578 


E 


3 


061110B 


21:35:40.400 


+06:52:33.91 


26.04+0.29 


> 21.24 


3.4344 


A* 


3 


061121 


09:48:54.564 


-13:11:43.09 


22.75+0.03 


20.14+0.19 


1.3145 


A* 


3 


070103 


23:30:13.793 


+26:52:35.32 


24.18+0.14 


21.07+0.36 


2.6208 


E 


13 


070110 


00:03:39.268 


-52:58:27.12 


25.02+0.11 


22.41+0.50 


2.3521 


AE* 


3 


070129 


02:28:00.915 


+11:41:04.56 


24.39+0.12 


21.01+0.18 


2.3384 


E 


13 


070224 


11:56:06.637 


-13:19:48.33 


25.96+0.31 


> 21.20 








070306 


09:52:23.310 


+10:28:55.28 


22.90+0.08 


19.70+0.10 


1.49594 


E 


39,16 


070318 


03:13:56.800 


-42:56:46.25 


24.43+0.11 


> 21.04 


0.8397 


A 


3,40 


070328 


04:20:27.615 


-34:04:00.63 


24.43+0.13 


> 20.98 








070330 






> 26.15 


> 21.19 








070419B 


21:02:49.801 


-31:15:48.82 


25.23+0.20 


> 21.22 


1.9591 


E 


13 


070506 


23:08:52.392 


+10:43:21.00 


26.10+0.22 


> 21.36 


2.3090 


A 


3 


070611 






> 27.09 


> 21.75 


2.0394 


A 


3 


070621 


21:35:10.061 


-24:49:02.18 


25.77+0.23 


> 21.51 








070721B 


02:12:32.950 


-02:11:40.80 


27.53+0.44 


> 22.38 


3.6298 


A* 


3 


070802 


02:27:35.722 


-55:31:38.76 


25.11+0.21 


21.27+0.41 


2.4541 


A* 


41,3 


070808 


00:27:03.310 


+01:10:35.81 


26.71+0.33 


21.77+0.37 








050412 


12:04:25.460 


-01:12:00.05 


25.55+0.19 


> 21.40 








050603 






> 26.61 


> 22.66 








060117 






> 26.28 










060223A 






> 26.25 


> 21.24 


4.406 


A 


28 


060313 








> 21.13 








060505 


22:07:03.500 


-27:48:55.57 


17.90+0.00 


16.22+0.01 


0.0889 


E 


42,43,44 


060602A 


09:58:16.750 


+00:18:12.56 


23.88+0.08 


20.56+0.27 


0.787 


E 


45 


070721 A 


00:12:39.048 


-28:33:00.14 


25.16+0.17 










070721A 


00:12:39.139 


-28:33:00.98 


23.14+0.15 










070810A 






> 26.73 


> 21.93 


2.17 


A 


46 


080207 


13:50:02.997 


+07:30:07.44 


25.84+0.18 




2.0858 


E 


13 



Note. — Columnl: GRB name; Columns 2 and 3: host coordinates; Column 4: R-band Vega magnitude; 
Column 5: ii's-band Vega magnitude; Column 6: redshift; numbers in boldface were determined from our 
program; Column 7: source of the redshift: (A)bsorption features in the afterglow spectrum (*, with fine- 
structure transitions), (E)mission lines from the host, or (P)hotomctric redshift of the afterglow; Column 
8: redshift reference. The horizontal line separates the bursts in the TOUGH sample from 10 extra hosts 
that were observed within our program. 
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Table 4. Median Properties of TOUGH and TOUGH Subsamples 



Sample 


R 


K s 


(2) low 


Median z 


( 2 >high 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


TOUGH 


25.52 ±0.23 


> 22.53 


1.50 




2.53 


FORS 


23.59 ±0.51 


20.78 ± 1.30 




1.25 ±0.29 




X-shooter 


24.48 ±0.15 


21.13 ±0.27 




2.24 ±0.15 




OA 


25.80 ±0.31 


> 22.53 


1.95 




2.32 


XRT 


24.36 ± 0.40 


21.13 ±0.29 








z 


25.11 ± 0.31 


> 22.53 




2.14 ±0.18 




No z 


25.88 ±0.24 


> 21.77 








R 


24.73 ±0.21 


22.15 ±0.69 


1.29 




2.39 


No R 






3.07 




3.97 


K s 


24.27 ± 0.28 


20.74 ±0.16 


1.10 




1.94 


No K a 


26.61 ± 0.43 




2.14 




3.35 



Note. — Column 1: TOUGH sample or subsample; Column 2: median 
R magnitude of subsample; Column 3: median K s magnitude of subsample; 
Column 4: minimum median rcdshift of subsample, assuming all unknown 
rcdshifts are less than the minimum redshift; Column 5: median redshift; 
Column 6: maximum median redshift of subsample, assuming all unknown 
rcdshifts are larger than the maximum redshift. 

Errors in the median are estimated as 1.4826 X MAD/^/N — 1, where MAD is 
the median absolute deviation, N is the number of systems used to determine 
the median and 1.4826 is a scale factor that applies to Gaussian distributions. 
The subsamples are the following: FORS: the six well-defined FORS emission- 
line redshifts, i.e., excluding the dropout rcdshifts (GRBs 050822, 051006, 
051117B, 060908, 061021, and 070306); X-shootcr: the nine new X-shootcr 
redshifts; OA: systems with an optical/NIR afterglow; XRT: systems with 
no optical/NIR afterglow; z: systems with a known redshift; No z: systems 
without a known rcdshift; R: systems with a measured R magnitude; No R: 
systems without a measured R magnitude (including the high-rcdshift systems 
GRB 050904 and GRB 060927 which were not observed); K a : systems with a 
measured K s magnitude; No K s : systems without a measured K s magnitude 
(including the high-redshift GRB 050904 but excluding GRB 050819 and GRB 
051006 which were not observed). 



